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AMENDMENTS TO THE CLAIMS: 

This listing of claims will replace all prior versions, and listings, of claims in the 
application: 

Claim 1 (Canceled). 

2. (Previously Presented) A method as in claim 4 wherein the fractal subsurface 
model is generated by a modified Voronoi tessellation technique which comprises modifying the 
Voronoi tessellation by using L p norm, where p is an exponent which can assume fractional real 
numbers and negative real numbers, thereby enabling greater range for generation of different 
fractal sub-surfaces, which are closer to a real geological situation. 

3. (Previously Presented) A method as in claim 4 wherein the natural setting of the 
geological subsurface being modeledjs selected from the group consisting of: a sedimentary 
basin, hydrocarbon deposits, oil reservoirs, aquifers and mineral deposits. 

4. (Currently Amended) A method for the generation of a three dimensional fractal 
subsurface structure by Voronoi tessellation and computation of gravity response of such fractal 
structure, by generation of fractal subsurface structures and computation of a forward gravity 
response of such structure for delineation and modeling of an underlying anomalous object, said 
method comprising: 

(a) selecting Voronoi centers at a plurality of locations over a region of interest, the 
Voronoi centers being represented by x, y, z co-ordinates; 

(b) generating an initial model of a subsurface fractal geological object, having variation 
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in a physical property in lateral and vertical directions; the initial model being generated by 
tessellating the Voronoi centers, and assigning values of physical property variations during 
generation of the model on the basis of pre-determined assumptions in the model; 

(c) assigning different discernable representations to regions in the model which have 
different physical properties to enable demarcation; 

wherein the computation of a gravity response due to a fractal subsurface generated by 
modified Voronoi tesselation is carried out by an analytical expression which comprises: 

(a) demarcating boundaries of regions having a different physical property in the 
tessellated region, the boundaries forming a polygonal shape in 2-dimensional space; 

V = Gp±[W arc cosfe / n )(x i+l I r M ) + {y t I n ){y i+l I r M )} 



i=l 



1/2 



-arc sin — = — + arc sin 

(A +^ ) (Pi +z t ) 

Where S= +1 if pi is positive, S= -1 if pi is negative, 

W= +1 if mi is positive, W= -1 if mj is negative, 

'Z' is depth and 'n' is number of sides in the polygon. 

G is universal gravitational constant, p is the density of the tessellated regions; 
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pi = perpendicular to a side of an irregular geometrical body from a point at which 
anomaly is being calculated; 

qi = cos(0O, where Q { is an angle between a side of the irregular body and line 
joining the first point (A) of the arm from the point at which anomaly is being computed; 

fi = cos((pi) 5 where (pi is an angle between a side of the irregular body and line 
joining the end point (B) of the arm from the point at which anomaly is being computed, 

mj = angle subtended by an arm of the irregular body at the point where anomaly 
is being computed; 

^ = distance of a point at the boundary of irregular body from the point at which 
anomaly is being computed; 

r i+1 = distance of next consecutive point at the boundary of the irregular body from 
the point at which anomaly is being calculated; 

r i>i+1 = length of the side of the irregular body; 

i = is a subscript, which is used to select the next consecutive point on the surface 
of the body; 

where the effect of a common arm of the adjacent polygon is removed; 
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(b) repeating the gravity response computation due to the polygon for all adjacent 
polygons of different physical properties using the demarcated polygon boundary; 

(c) adopting the process of gravity response computation for tessellated regions lying at 
different depths to obtain a response; and 

(d) integrating the response using Simpson/Gauss quadrature formula at plurality of grid 
nodes overlain on the region of interest ; and 

(e) creating a model of the underlying anomalous object . 

5. (Previously Presented) A method as claim 4 wherein results are achieved through 
self-written software, which generates a Voronoi tessellated subsurface region and computes a 
gravity response of the same. 

6. (Previously Presented) A method as claim 4 which is applied to geophysical 
inversion, wherein the tessellated regions are altered by changing the position of Voronoi centers 
in each iteration. 
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